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The detailed formulation for loop quantum cosmology (LQC) in the Bianchi I model with a scalar 
massless field has been constructed. In this paper, its effective dynamics is studied in two improved 
strategies for implementing the LQC discreteness corrections. Both schemes show that the big bang 
is replaced by the big bounces, which take place up to three times, once in each diagonal direction, 
when the area or volume scale factor approaches the critical values in the Planck regime measured 
by the reference of the scalar field momentum. These two strategies give different evolutions: In 
one scheme, the effective dynamics is independent of the choice of the finite sized cell prescribed 
to make Hamiltonian finite; in the other, the effective dynamics reacts to the macroscopic scales 
introduced by the boundary conditions. Both schemes reveal interesting symmetries of scaling, 
which are reminiscent of the relational interpretation of quantum mechanics and also suggest that 
the fundamental spatial scale (area gap) may give rise to a temporal scale. 
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I. INTRODUCTION 

The comprehensive formulation for loop quantum cos- 
mology (LQC, see [T] for a review) in the spatially flat 
and isotropic model has been constructed in detail [2 [3], 
giving a conceptual framework and a solid foundation to 
analyze the physical issues of the quantum theory. With 
a massless scalar field serving as the emergent time, the 
investigation shows that the quantum evolution is deter- 
ministic across the deep Planck regime and in the back- 
ward evolution of the states which are semiclassical at 
late times, the big bang is replaced by a big bounce. 

In the original construction (called "/i -scheme" ) of 
[2J [3J , the discreteness variable introduced to impose the 
fundamental discreteness of quantum geometry is taken 
to be constant (referred to as fi ). However, it has been 
shown that this prescription leads to the wrong semi- 
classical behavior that, in particular, the critical value of 
matter density at which the big bounce takes place can 
be made arbitrarily small by increasing the scalar field 
momentum p^. To fix this problem, the construction 
was further improved by a more sophisticated implemen- 
tation of the underlying physical ideas of loop quantum 
gravity (LQG) with the discreteness variable (referred to 
as /2) varying adaptively [41] . In this improved dynam- 
ics (called "/i-scheme" ) , the big bounce occurs precisely 
when the matter density enters the Planck regime, re- 
gardless of the value of p^. This construction was also 
extended to k = +1 and k = — 1 Robertson- Walker mod- 
els [SJE]. 

In order to further develop this formulation and ex- 
tend its domain of validity, based on the same principles 
of [2J OS], both the precursor strategy (^-scheme) and 
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the improved strategy (/S-scheme) were applied and re- 
constructed for the Bianchi I model to include anisotropy 
[TJ. The analytical investigation shows that the state in 
the kinematical Hilbert space associated with the clas- 
sical singularity is completely decoupled in the difference 
evolution equation, indicating that the classical singu- 
larity is resolved in the quantum evolution and the big 
bounce may take place when any of the area scales un- 
dergoes the vanishing behavior. 

While a thorough numerical investigation for [7J re- 
mains to be done to draw the definite conclusion for the 
details of the quantum evolution in the Bianchi I model, 
the effective dynamics with LQC discreteness corrections 
was first studied in [5] using ^-scheme for the vacuum so- 
lution and the /i-scheme effective dynamics has also been 
done in [9] for the case with a massless scalar field and 
in [TU] with the inclusion of generic matters. Not only 
do the results at the level of effective dynamics agree 
with the anticipations in [7] but more intuitive pictures 
are also obtained in the semiclassical approach, giving 
insights into how and why the big bounces take place. 

Because of the variety of Bianchi I models, extend- 
ing the improved scheme is slightly ambiguous as a few 
strategies are possible. Among them, two schemes are 
of particular interest (referred to as "/1-scheme" and u pf- 
scheme" in this paper). The /i-scheme is suggested in [7J 
and used in [51 [TU] , while the r/-scheme is disfavored by 
[7] and only mentioned briefly in the appendix of [TU] . It 
is later realized that the virtues of //'-scheme might have 
been overlooked and we should not dismiss either of them 
prematurely in favor of the other before more physics is 
understood, as either scheme has both advantageous and 
disadvantageous features. 

To better understand the difference of their phe- 
nomcnological ramifications, this paper presents the de- 
tails of effective dynamics for the Bianchi I model with 
a massless scalar field in both ft- and /i'-schemes. The 
investigation shows that in both schemes the big bang 
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singularity is replaced by the big bounces, which take 
place up to three times, once in each diagonal direction. 
However, p- and /2'-schemes give different evolutions, dis- 
tinct from each other not only in detail but also qualita- 
tively. For instance, the indication of the occurrence of 
big bounces is the matter density in /2'-scheme while 
it is the "directional density" qj in /1-scheme. Further- 
more, in the different schemes, the big bounces conjoin 
the different pairs ("antipodal" pair vs "conjugate" pair 
) of classical solutions on the two asymptotic sides of 
classical regime. 

In the Hamiltonian framework, we have to restrict the 
spatial integration to a finite sized cell V to make the 
Hamiltonian finite. This prescription raises the ques- 
tion whether the resulting dynamics is independent of 
the choice of V. This issue is very subtle and has caused 
a great deal of confusion because the dependence on V 
of the Ashtekar variables is rather obscure. In order to 
elucidate the subtleties, the classical dynamics is care- 
fully studied both in metric variables and in Ashtekar 
variables; the precise physical meanings of Ashtekar vari- 
ables are then obtained. With the subtleties clarified, it 
can be shown that the effective dynamics in /i'-scheme is 
completely independent of the choice of V as is the clas- 
sical dynamics, while the effective dynamics in /1-scheme 
in fact reacts to the macroscopic scales introduced by the 
boundary conditions of V. This is an important differ- 
ence between these two schemes. 

In addition to the issues related to the dependence on 
V, the effective dynamics in both schemes also reveals 
other interesting symmetries of scaling, which are rem- 
iniscent of the relational interpretation of quantum me- 
chanics. It is also suggested that the fundamental scale 
(area gap) imposed for the spatial geometry may gives 
rise to a fundamental scale in temporal measurement. 
These observations could simply be technical artifacts 
but encourage further investigations for more sophisti- 
cated models to tell their profundity. 

This paper is organized as follows. In Section [H] the 
classical dynamics of the Bianchi I cosmology with a 
massless scalar source is solved in Hamiltonian formal- 
ism in terms of both metric variables and Ashtekar vari- 
ables. The effective dynamics with LQC discreteness cor- 
rections is constructed and solved in Section IIIII for the 
p- and /l'-schemes respectively. The scaling symmetry 
and issues about relational measurements are discussed 
in Section |IV| Finally, the results and outlooks are sum- 
marized in Section [VJ As a comparison to p/ p'-scheme, 
the effective dynamics in /i -scheme is also included in 
Appendix |Aj The details of heuristic arguments and mo- 
tivations for p- and /^'-schemes are given in Appendix |B| 



II. CLASSICAL DYNAMICS 

Before considering the LQC discreteness corrections, 
in this section, we study the classical dynamics of the 
Bianchi I cosmology with a massless scalar field. In order 



to grasp the precise physical meanings of the Ashtekar 
variables, which are the fundamental degrees of freedom 
in LQC, the formulations in metric variables (also see 
[TT]) and in Ashtekar variables are both investigated. 



A. Formulation in metric variables 

The Hilbert-Einstcin action for the cosmology with a 
massless scalar field <f> minimally coupled to the gravity 
is given by 



S = 5 grav + Scf, := I d 4 x £ := I dt L 



(2.1) 



iLS dix ^ R+ \! d 



In Bianchi I models, the spacetime metric can be written 
in the diagonal form: 



ds 2 = g^dx^dx" 



(2.2) 



-dt 2 + a\{t)dx 2 + al(t)dy 2 + aj(t)dz 



and the homogeneous field cj>(x, t) = cf>(i) is independent 
of the spatial coordinates. The temporal variable t is 
the proper time and aj are the diagonal scale factors. 
The coordinates (x,y,z) which diagonalize <? M „ are the 
co-moving coordinates. 



The metric (2.2) gives the Ricci scalar: 



R = 2 ^l + ^ + ^ + ^l + ^ + ^l). (2 .3) 
\a 1 a 2 a 3 a 2 a 3 a x a^ a x a 2 ) 

This leads to 

Sgrav = y d 4 x|-^ (dia 2 a 3 + aia 2 a 3 + aia 2 d 3 ) 

-aid 2 d 3 - dia 2 d 3 - di^fflaj. (2.4) 

The total time derivative can be ignored and we then 
have the Lagrangian 

L = I d z x\ -— *— (aid 2 d 3 + dia 2 d 3 + dxd 2 a 3 ) 

57TG 



-aia 2 a 3 - 



(2.5) 



Because of homogeneity, the spatial integration diverges 
as we are considering the noncompact Bianchi I model. 
To circumvent this problem, the integration is restricted 
to a finite sized cell V; that is, we take 



d 3 x = dx I dy dz =: V, (2.6) 
'v Jo Jo 

where Lj are the coordinate lengths of the edges of V 
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and V is its coordinate volume. 1 [This prescription is 
equivalent to compactify the space or impose a spatial 
periodicity. We will see that this does not change the 
classical dynamics but might have effect on the quantum 
corrections.] 

With this prescription, the Lagrangian is given by 



L = V<- 



1 



8nG 



(aid 2 d 3 + 0,10-20,3 + did 2 a 3 ) 



+aia 2 a 3 - 



(2.7) 



and the canonical momenta of <p and a/ are defined as: 

(2.8) 



■ = — r = Vaia 2 a 3 <p, 



7Tl 



dL 
dai 



V 



(d 2 a 3 + a 2 d 3 ) 



(2.9) 



and so on for 7r 2 , 7r 3 in the cyclic manner. 2 [Note that the 
definitions of momenta depend on the cell V we choose!] 
Solving (2.9 1 gives the velocities: 



«i 



AllG ( 0\~K\ 7T 2 7T 3 



V \a 2 a 3 
by which we find the Hamiltonian 



0-2 



H 



TTiai + 7r 2 a 2 + 7r 2 a 2 



(2.10) 



(2.11) 



V~ 



2vrG 



ai7r 2 



27T 2 7r 3 27Tl7r 3 27Tl7T 2 



a i 



a 2 



«3 



a 2 7r| a^nl 1 



a 2 a 3 aia 3 
with the canonical relations: 



a\a 2 



2aia 2 a 3 



{0,P4>} = !) 
{a/,7rj} = (5/,/. 



(2.12) 
(2.13) 



We can simplify the Hamiltonian by choosing the lapse 
function ./V = Va\a 2 a 3 and thus introducing the new 
time variable dt 1 = (Vaia 2 a 3 )~ 1 dt. The rescaled Hamil- 
tonian is then given by 

H = 2irG{— 2a 2 a 3 7r 2 7r 3 — 2aia 3 7Ti7r 3 — 2aia 2 7Ti7r 2 



, 2 2 , 2 2 , 2 2\ , 

+a 1 TT 1 + a 2 7r 2 + a 3 7r 3 ] + 



(2.14) 



1 The coordinate lengths Lj and volume V are independent of time 
once the cell V is chosen since we are working on the co-moving 
coordinates. The physical lengths of the edges and volume of V 
are given by Lj := ajLj and V := a\a2a^V respectively, which 
evolve with time. 

2 We will not mention the obvious cyclic repetition hereafter. 



The equations of motion are governed by the Hamilton's 
equations: 



dp^ 
~dV 
d(j> 
dt 1 
dai 

w 

dni 
~dY 



{<f>,H'} = 
{ax,H'} 







dH' 



is constant, (2.15) 
(2.16) 



= AirGai {ai^i — a 2 ir 2 — a 3 7r 3 ) , 



dai 



= -A-kG-ki (a 1 ir 1 - a 2 n 2 - a 3 7r 3 ) . 
In addition, the Hamiltonian must vanish: 
if'(a/, 77/,^) = 0. 



Combining (2.171 and (2.181 gives 
d 



7 (a/7r/) = 0. 



We assign 



dt 



o\t\\ = -h{K 2 + /C 3 ) 



(2.17) 
(2.18) 
(2.19) 

(2.20) 
(2.21) 



with the dimensionless constants /C/, which will be used 
to parameterize the solutions of evolution. Equations 
(12.16} and (12.17} then give 



1 dai 

ai d(f> 



InGh 



V&^G 



if we define 



:= hV&KGK, 



(2.22) 



(2.23) 



with the dimensionless constant K,^. Regarding <f> as the 
emergent time, the solutions of evolution are given by 



a/(</>) = aj(</>o) e 



(2.24) 



where we scale the parameters ICi — ICkj, JC<p = ICn^ 
(with JC > 0) such that 



Kx + K 2 + K 3 = ±1. 



The Hamiltonian constraint (2.19) now reads as 



= 2 (/C 2 /C 3 + /Ci/C 3 + K,\K, 2 ) 

or equivalently 



1. 



(2.25) 



(2.26) 



(2.27) 



The classical Bianchi I model with a massless 
scalar field admits both "Kasner-like" (two of kj 
positive/negtive and the other negative/positive) and 
"Kasner-unlike" (all K/ positive/negative) solutions. The 
Kasner-like solution, which has two expanding and one 
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contracting directions (say, > and K1+K2 + K3 = 1), 
eventually encounters the "Kasner-like singularity" (a 
given regular cubical cell stretches as an infinitely long 
line) in the far past and the "planar collapse" (a regu- 
lar cubical cell stretches as an infinitely large plane) in 
the far future. On the other hand, the Kasner-unlikc 
solution, with all directions expanding, encounters the 
"Kasner-unlike singularity" (a regular cubical cell van- 
ishes to a point) in the far past and no planar collapse. 3 

We will see that with LQC discreteness corrections, 
both Kasner-like and Kasner-unlike singularities are re- 
solved and replaced by the big bounces, whereas the pla- 
nar collapse remains its destiny even one of the three di- 
agonal directions approaches infinitely small length scale. 

Note that, by (2.8 1 and (2.9), and 717 depend on the 



choice of V and scale as p^, 717 oc V. Thus, the directly 
measurable quantities are not the canonical momenta but 
rather the momentum density of 0: 



where, again, the dependence of V is gone. Therefore, 
whether in terms of the proper time t or the emergent 
time 0, the classical dynamics is completely independent 
of the finite sized cell V we choose to make sense of the 
Hamiltonian formalism. [The independence of the choice 
of V is not necessarily retained when quantum correc- 
tions are taken into account. In any case, however, the 
dynamics is independent of V in terms of t if and only if 



it is so in terms of 0.] Furthermore, (2.33 1 gives a simple 
relation: 



t-to = 



K<j> 



e * — 1 



(2.34) 



for the classical solutions given by (2.241, but the rela- 
tion between t and is more complicated if the quantum 
corrections are taken into account. 



^ = = V _1 W87rG£*, (2.30) 



B. Formulation in Ashtekar variables 



which is nothing but the time derivative of 0, and the 
quantities: 

= " 27g {H2 + Hs) = - V (/C2 + Kz) (2 ' 31) 

with the directional Hubble rates defined as 



H, 



(>i 



Li' 



(2.32) 



For given initial physical conditions 0| to and Hj\ toJ the 
constants of motion K.^ and /Cj both scale as cx V| to , the 
physical volume of V at the initial time to. The ratio 
ICi/K-cf, — Ki/Kc/, is nevertheless independent of V and 
hence the solution of ai(<j>) given by (2.241 does not de- 
pend on the choice of V. Once a/(0) is solved, to know 
the solution aj(i) as a function of t, we only need to 
convert back to t via 



t-tn = 



1 dt= t* Va^a^g^) 

t J<t>o P<t> 

ai (0)02(^)03(0) d , 
4>o ai(0o)fl2(0o)a3(0o) 0| to 



(2.33) 



3 In the case of no matter sources, the classical dynamics of vacuum 
Bianchi I model yields the standard Kasner solution: 



ai(t) = ai (to) ( 



where ki,K2,^3 are constants subject to: 

rei + + «3 = 1 and k\ + k| + k 



(2.28) 



(2.29) 



which are to be compared with l |2.25| l and l |2.27[ l. Apart from 
the trivial solution (where two of kj vanish), all Kasner solutions 
must have two of kj positive and the other negative. 



Equivalently, the classical dynamics studied above can 
be reformulated in terms of Ashtekar variables as the 
groundwork for LQC formulation. We follow the lines 
of Section III of [7 but emphasize the explicit relations 
between the metric and Ashtekar variables by rigorously 
taking care of the subtleties overlooked (and correcting 
the mistakes made) in [71 [3] ■ 

In Bianchi (homogeneous) models of cosmology, on the 
homogeneous space-like slice S, we can choose a fiducial 
triad of vectors °ef and a fiducial co-triad of covectors °oj l a 
that are left invariant by the action of the Killing fields 
of E. (Note °ef%l = 8%.) The fiducial 3-metric of E is 
given by the co-triad °uj l a : 



°<Z„ b = <°^%. (2.35) 



In LQG, the Ashtekar variables consist of the canon- 
ical pairs: the densitized triads E i(x) and connections 
A a l (x), which satisfy the canonical relation: 

{A a \x),E b 3 (x')} = 8wG 7 6 i j S b a 6 3 (x-x'), (2.36) 

where 7 is the Barbero-Immirzi parameter. In the con- 
text of Bianchi models, it is more convenient to consider 
the reduced canonical variables: 



AJ := 



d 3 xV^(°q- 1/2 °to l a E a 3 ) 



d 3 x% l a E a 3 , 



d 3 xV^{Aj°et) 



(2.37) 
(2.38) 



The reduced canonical variables are essentially the same 
as the original canonical variables, but stripped of the 
spatial dependence enforced by the Bianchi symmetry 
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and integrated over a finite sized cell V. 4 Conversely, we 
have the replacement rules: 



V-'V^E^ °e? = V~ X EU °el (2.39) 



T/ _1 A °rJ 



(2.40) 



where V is the fiducial volume defined as the volume of 
V measured as if the 3-metric was °q a b, i.e. 



d 3 Xy/^. 



The canonical relation (2.361 now reads as 



{A i j ,E k l } = 8nG 1 V 6%Sj. 



(2.41) 



(2.42) 



^,6 



Howe ver, if we com pute {A a l ,E j} using (2.39) and 
(2.401 through (2.42), the 5-function in the right hand 



side of (2.361 is correspondingly replaced with 
6 3 (x-x') — > V- 1 



(2.43) 



as expected due to the regularization procedure. 

For the diagonal Bianchi class A models (Bianchi I 
included), E l j and A? are diagonalizable and their di- 
agonalized entries are denoted as pi and 5/ (J = 1,2,3) 
respectively. Equation (2.421 then becomes 



(2.44) 



In Bianchi I models, correspondingly, we choose the fidu- 
cial 3-metric to be the diagonal form: 



o„2 





°q ab = \ °a\ 
°a\ 



(2.45) 



and the finite sized cell V is adapted to be rectangu- 
lar with respect to °g a b and coordinated as described in 
dZ6|. The fiducial volume of V is now V Q = Li°L 2 °L 3 = 



ai°a2°azV with 

°L/ := / ° ai d Xl = °a I L I (2.46) 

JO 

being the fiducial lengths of the edges of V. 5 



4 The spatial integration over a finite sized cell V is not explic- 
itly noted in [7] and nor in most literatures. This prescription 
is subtle but important because otherwise the Poisson bracket 
{-E 1 j , Ai J } would be distributional. 

5 The fiducial length/ volume should not be confused with the co- 
ordinate or physical length/volume (see FootnoteJTJ. (Also note: 
Vo and °Lj could depend on t, depending on the choice of °q a b-) 
Most literatures do not distinguish these notions completely, but 
the nuances are highlighted in this paper in order to give the 
precise physical interpretations of cj and pj. 



To get rid of the nonphysical dependence on °q a b, it is 
convenient to introduce the new variables: 



Pi ■= "L^pi, 

ci := CWLi)- 1 ^, 

which satisfy the canonical relation: 

{ci,Pj} = 8ttGjSij. 



(2.47) 
(2.48) 



(2.49) 



The relation between the densitized triad and the 3- 
mctric is given by 



~ a ~b 



qq ab = S l 3E iE j, 



which leads to 



Pi = L 2 L 3 (a 2 a3) = L 2 L 3 



(2.50) 



(2.51) 



Thus, the triad variables pi are the physical areas of the 
rectangular surfaces of V. 6 Comparing the canonical re- 
lations (12 131) and (12. 491) via (12. 5 1L we have 



4irjG 
L 2 L3a2a3 



(a 1 n 1 - a 2 7r 2 - a 3 7r 3 ) 



By (2.9), we have 



(2.52) 



(2.53) 



which tells that, classically, the connection variables c/ 
are the time change rates of the physical lengths of the 



edges of V (up to the constant 7). Note that (2.51 1 and 



(2.521 remain the same even with the quantum correc- 
tions since the symplectic structures given by fl2T3) and 
(2.491 as well as the relation (2.50) arc unchanged. By 



contrast, (2.531 is modified in the presence of quantum 



effects because Hamiltonian is changed and thus (2.9) no 
longer holds. 



By (2.511 and (2.521, the classical Hamiltonian (2.111 



can be recast in terms of c/ and pj: 



H = H 

(c 2 P2C 3 p 3 



H<t> 

- cipxc 3 p 3 



(2.54) 



C1P1C2P2) 



Pi 



inG^ 2 yjpxpipz 



2^PlP2P3 



[We can also derive i? gra v directly from the Hamiltonian 
constraint of the full theory of LQG. For this approach, 
see the text toward (B4| in Appendix [B] or Section III.C 
of [7].] The solutions of cj and pi can be readily ob- 



tained by translating (2.211 and (2.24) via (2.511 and 
(2.521. But in order to draw a parallel to solving the ef- 



fective dynamics with LQC corrections studied later, in 
the following, we present solving the equations of motions 
directly for cj and pj . 



6 More precisely, | |2.5l| should be |pi| = L2L3. The sign of pj 
corresponds to spatial orientation, which we ignore in this paper. 
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Again, we simplify the Hamiltonian by choosing the 
lapse function N — ^/p\P2P3 — Vaia,2a 3 and thus in- 
troducing the new time variable dt' = {p\P2P3) 1 ^ 2 dt = 
(^010203) dt. The rcscalcd Hamiltonian (which is the 
same as (2.14l) is given by 



H' = - 



(c 2 p 2 C 3 P3 + C1P1C3P3 + C1P1C2P2) 

8nGj 2 



(2.55) 



The equations of motion are governed by (2.15) and 
(2.161 as well as 



dc\ 
~dt> 

dpi 
dt' 



{c l ,H'} = 8nG 1 d ^- 
opi 

-7 _1 ci (c 2 P2 + c 3 p 3 ) , 
{p l ,H'} = -^G 1 d]l 



dci 



7 X pi (c 2 p 2 + C3P3) • 



(2.56) 



(2.57) 



In addition, the constraint that the Hamiltonian must 
vanish yields 



ff'(c/,p/,p0) = 



(2.58) 



(C2P2C3P3 + CiP!C 3 p 3 + C1P1C2P2) 



[Note that substituting ( |2.51[ ) into ( |2.57[ ) leads to ( |2.53[ ) 
again.] 



Combining (2.561 and (2.57) gives 
d 



dt 



-,{PiCi) = 0, 



PtC! = 8nG-fhlC I , (2.59) 



where ICi are the same constants as defined in (2.211 



Taking (2.591 into (2.581 yields the same constraint 
(|2~26l. 



Putting (2.59) into (2.571 gives 



1 dpi 
p 1 dt' 



inGh(IC2+IC 3 ). 



Referring to (2.161, this leads to 



Pi d(p p^ \ 

and consequently 

PiW =Pi{<Po)e \ K 4> J 



(2.60) 



(2.61) 



(2.62) 



which is the same solution as given by (2.241 



III. EFFECTIVE DYNAMICS WITH LQC 
DISCRETENESS CORRECTIONS 



In LQC, the connection variables c/ do not exist and 
should be replaced with holonomies. In the effective the- 
ory, to capture the quantum corrections, following the 



procedures used in the isotropic case 
prescription to replace cj with 



cr 



sin(/i/c/) 



/./,/ 



J, we take the 



(3.1) 



introducing the discreteness variables pi. The heuristic 
argument starting from the full theory of LQG and lead- 
ing to (3.1 1 can be found in Appendix |B) This prescrip- 



tion can also be understood as the WKB approximation 
of the full quantum theory of LQC [8 . 

Imposing this prescription plus the loop quantum cor- 
rection to the inverse triad on ( 2.54[ ) , we have the effective 
Hamiltonian constraint to the leading order: 



H cS = /(pi)/(pa)/(P3) 



V% f(Pl)f(P2)f( P 3) 



(3.2) 

P2P3 + cyclic terms 



2 8ttG7 2 
sin(^ 2 c 2 ) sin(/2 3 c 3 ) 

p2p-3 

where f(pi) is the eigenvalue of the inverse triad oper- 
ator 1 / y/pj. The loop quantization gives the quantum 
corrections on f(pi): 



f(Pi) 



oc p/"/4? +1 for PI « 4, ( ' 



with the Planck length £p\ := v Gh and a positive n. 
The corrections to f(pi) are significant only in the deep 
Planck regime in the vicinity of pi — 0. From now on, 
we will ignore the quantum corrections to f(pi) by sim- 
ply taking its classical function f(pi) = pi^ 1 ^ 2 . The fact 
that this effective quantum modification is a good ap- 
proximation for the sates which are semiclassical at late 
times has been verified in the isotropic models of LQC 
sourced with a massless scalar field [H [13]. The valid- 
ity of this effective prescription for the vacuum Banchi 
I model was discussed in the context of WKB approx- 
imation in [S] and that for the Bianchi I model with a 
massless scalar field can also be affirmed [14 j. [We will 
see that in the backward evolution of the solutions which 
are semiclassical in the late times, the big bounces take 
place much earlier before the discreteness correction on 
the inverse triad operator becomes considerable, and it is 
the "non- locality" effect (i.e., using the holonomies) that 
accounts for the occurrence of big bounces.] 



With f( Pl ) 



Pi 



-1/2 



by choosing N = y/p a P2P3 



and dt' = (P1P2P3) 1 ^ 2 dt, the Hamiltonian (3.2 1 can be 
rescaled as 



(3.4) 



8ttG 7 2 



sin(/i2C 2 )sin(/x 3 C3) 



P2P3 + cyclic terms 



As for imposing the fundamental discreteness of LQG 
on the formulation of LQC, the precursor construction 
(/i -scheme) is to take /2/ as constants (referred to as [if). 
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However, it has been shown that in the isotropic case 
/i -scheme can leads to the wrong semiclassical limit and 
should be improved by a more sophisticated construction 
(/2-scheme) in which the value of p, depends adaptively 
on p (via M oc 1/^Jp) and thus implements the underly- 
ing physics of quantum geometry of LQG more directly 
[1]. Extending this scheme to the Bianchi I model is 
slightly ambiguous as a few possibilities exist. Among 
them, two constructions are most preferable (referred to 
as "/2-scheme" and "M'-scheme" ) : 



/i-scheme: 




M3 



/i'-scheme: 



/ PiA 
P2P3 



M 2 



/ P2A 
PWa 



M 3 



/ p 3 A 
P1P2 



(3.5) 



(3.6) 



Here A = ^(4777^) is the area gap in the full theory 
of LQG. 

Either scheme has its own advantages and drawbacks 
and until more detailed physics is investigated it is still 
debatable which one makes better sense. In particular, 
M-scheme is suggested in [7], since in the construction 
for the full theory of LQC the Hamiltonian constraint in 
/2-scheme gives a difference equation in terms of affinc 
variables and therefore the well-developed framework of 
the spatially flat-isotropic LQC can be straightforwardly 
adopted. (However, it is argued in [15 that /2-scheme 
may lead to an unstable difference equation.) By con- 
trast, M'-scheme does not admit the required affine vari- 
ables and the full LQC of it is very difficult to construct. 
On the other hand, /2'-scheme has the virtue that its ef- 
fective dynamics is independent of the choice of V as will 
be seen. (More detailed comparisons and the heuristic ar- 
guments for both schemes are presented in Appendix |B|) 
To explore their differences, we study both /2-scheme and 
M'-scheme in the context of effective dynamics in Sec- 
tion |III A| and Section III B respectively. (As a compari- 
son, the effective dynamics in /i -scheme is also presented 
in Appendix [X] where we sec the insensible behavior that 
pi can be made arbitrarily small.) 



~dF 



d H- 

{c 1 ,H' fl } = 8nG^-^ 

3sin(/2iCi) ciCOs(/2iCi) 
2fh 2 

sin(/2 2 c 2 ) sin(/2 3 c 3 ) 

z P2 H — 

M2 M3 



-P:i 



(3.9) 



dp 



} = { Pll H' p } = -8nG 1 dh; ' 



7 1 piCOs(/iiCi) 

"sin(/2 2 c 2 ) 



/'2 



as well as 



-P2 



Pi 



sin(/2 3 c 3 ) 
M3 



-Pa 



(3.10) 



(3.11) 



AttGj 2 



sin(/i 2 c 2 )sin(/i 3 c 3 ) 



M2M3 



P2P3 + cyclic terms 

0, we have 



[Note that in the classical limit /i/C/ 
sin(/2/c/)//2/ — * cj, cos(/2/c/) — * 1 and therefore (|3.9[ ) 



(3.111 reduce to their classical counterparts (2.56) 



(2.581.] Also notice that substituting (2.511 into (3.101 



gives 



sin(/i 1 ci) 
Mi 

7L1 f Hi + H 3 | Hi + H 2 ff; 

2 \cos(/2 2 c 2 ) cos(/2 3 c 3 ) cos(/2iCi) 



(3.12) 



which is the modification of (2.53) with quantum correc- 
tions. 



Combining (3.9 1 and (3.10), we have 



3sin(/i/cr) cj cos(/i/Cf) \ dpi ,_ .dci 

1 ' p/cos(/i/c/) — 



2/2/ 



dt' 



_d_ 



Pi- 



si n (/ijcj) 
Mi 



0. 



dt' 
(3.13) 



which, in accordance with the classical counterpart 
(E5l, yields 



Pi- 



sin(/i/c/) 



Hi 



8irGjhlCi 



(3.14) 



A. Effective dynamics in /^-scheme 

The effective dynamics in M-scheme is specified by the 
Hamiltonian (3.4 1 with M/ given by (3.5). Again, the 



equations of motion are governed by the Hamilton's equa- 
tions and the constraint that the Hamiltonian must van- 
ish: 



dt' 

dd> 
df 



p<jj is constant, (3.7) 
(3.8) 



Taking (3.141 into (3.11 ) again gives the same constraints 



on the constant parameters as in (2.261 



Substituting ( |3.14[ ) into ( |3.10[ ) yields 
1 dpi 



pi dt' 



8-TrGh cos(/iiCi)(/C 2 + K.3). 



(3.15) 



By regarding </> as the emergent time via (3.8) and ex- 



pressing cos a; = ±vl — sin 2 cc, (3.151 then leads to 

1/2 



1 dpi 

Pi dej) 



±V8^G 



1 — Ki 



Qi 



Ql crit 



(3.16) 
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where we define the directional density for the I- 
direction: 



by (2.301 and (2.511. Therefore, contrast to the classical 



Ql 



Pi" 



(3.17) 



and its critical value is given by the Planckian density 
Ppi times a numerical factor k?/ ' Kj. 



Qi crit 



p Ph pn-^iSirG^A)- 1 . (3.18) 



[Note that the directional density gi is of the same di- 
mension as the matter density p^ :— v\I^V\V2V'i) an d 
thus the name.l 



From (3.161, we conclude: The big bang singularity 



(whether Kasner-like or Kasner-unlike) is replaced by the 
big bounces, which take place up to three times, once in 
each direction whenever each of gi approaches its criti- 
cal value g/, cr it (and thus cos(/Z/c/) flips sign in ( |3.15 1). 
Furthermore, the differential equations of individual pi 
are completely decoupled from one another and as a re- 
sult the epochs of the big bounces in different directions 
can be arbitrarily separate, depending on the given ini- 
tial condition. Also notice that the constants kj remain 
the same before and after the bounces. Equivalcntly, it 
can be rephrased that the big bounces connect the classi- 
cal solution specified by m, k.2, ^3 on one asymptotic side 
with its "antipodal" counterpart on the other asymptotic 
side specified by the "antipodal" parameters: 



«1,«2)«3 



—Ki, —K 2 , —K3. 



(3.19) 



On the other hand, however, the planar collapse in the 
Kasner-like case is not resolved but the vanishing behav- 
ior of one of the length scale factors ai continues. 

For given initial conditions, the differential equation 
(3.161 can be solved numerically. 7 The evolutions of 
Pi{4>), aj((/)) and Qi{4>) are depicted in parts (a), (b) and 
(c) respectively in FIG. [I] for Kasner-unlike solution and 
in FIG.[2]for Kasner-like solution. Note that the bounces 
occur at the moments exactly when gi approach their 
critical values. Also notice that each individual curve of 
Pi(4>) and ai((f>) in logarithmic scale becomes a straight 
line before and after the bounce, with the same slope 
but opposite sign, indicating that the bounces conjoin 
the two antipodal classical solutions. 

It is noteworthy that the directional density gi is the 
indication of the bounces but the quantity of gj is not 
independent of the choice of V, as we know 



V 



2A2 



Pi 



^2^3 



T 2 

— — — ( 
L2L3 



(3.20) 



7 Equivalently, we can directly solve the coupled differential equa- 
tions ( |3 - 8| ) — ( f3 - 10 [ ) for given pi((f>o), ci(<j>o) and p^. Also note 
that, to give the initial condition, an alternative way is to spec- 
ify pi(4>o) together with the consta nts Ki and K^. Given p/(0o) 
and fCi, ci(4>o) are fixed by (3.14l. 



dynamics, the effective dynamics in /2-scheme does de- 
pend on the choice of the finite sized cell V. Another sub- 
tler dependence on V comes from the fact that for given 
initial physical conditions <p\t and Hi\ ta , the constant of 
motion /C^ scale as oc V| to but tCi scale as oc V| t(J only 
approximately (see (3.121 and ( |3.14 1 and notice (3.12) in- 
volves the quantum modification with terms cos(/i/C/)). 
As a result, the ratio IC^/JCi = k^/ki is slightly depen- 
dent on V. (Nevertheless, in the classical regime, we still 
have k^/ki « (SttG) 1 / 2 ^/.) 

The problem not to be invariant under different choice 
of V is absent in the //-scheme as will be seen in Sec- 
tion III B However, we should not dismiss /2-scheme im- 



mediately as it is a common phenomenon that a quan- 
tum system reacts to macroscopic scales introduced by 
boundary conditions (for instance, the well-known "con- 
formal anomaly" as a "soft" breaking of conformal sym- 
metry) . If we have good physical input to tell what ex- 
actly the space is to be enclosed as V (such as in the 
compactified Bianchi I model, or applied for the finite 
sized homogeneous patches in the Belinsky-Khalatnikov- 
Lifshitz (BKL) scenario), the dependence on V could be 
rather meritorious than problematic and the bounce oc- 
currence condition (gi — g/, cr it) can be understood as: 
The physical areas (p\ — L2L3, etc) of the surfaces of 
V get bounced when each of them undergoes the Planck 
regime (times a numerical value Kj/k^) measured by the 
reference of the momentum ru. 



B. Effective dynamics in //-scheme 

The effective dynamics in /i-scheme is specified by the 
Hamiltonian (3.4 1 with /2j replaced by p\ given in (3.6 1. 
For simplicity, we choose a different lapse function N — 
(P1P2P3) 1 an d correspondingly introduce the new time 
variable dt" = {p\P2P3) 1 ^ 2 dt. With the new lapse, the 
Hamiltonian (|3.4|) is further rescaled to the simpler form: 



H'l, = (3.21) 
P<f, - 87rG ^ 2A {sin(/2 2 c 2 ) sin(/i3C 3 ) + cyclic terms} , 



where p^ is the matter density of <fi: 



2P1P2P3 2p 3 



(3.22) 



and we define p 3 :— P1P2P3 for convenience. Because 
|sin(/ljC/)| < 1, the vanishing of the Hamiltonian con- 
straint H'L = immediately implies 



< 



8nGj 2 A 



3p P1 . 



(3.23) 



The fact that p^ is bounded above implies that the big- 
bang singularity is resolved and the big bounce is ex- 
pected to occur when the matter density approaches 
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Planckian density. The inequality of (3.23) is saturated 
for the perfectly isotropic universe (i.e. K\ — k 2 — K3 = 
1/3) but we shall see that the maximal value of is 
in general lower than 3/?pi with anisotropics introduced. 
(Also sec Appendix C of [TO]-) 

To know the detailed dynamics for each individual pi, 
in addition to the Hamiltonian constraint, we study the 
Hamilton's equations: 



dt" 

dt" 



{p^H"} = Q 



is constant, 



{<l>, Hp} = p<f,(piP2Pa) \ 
dH" 



= {c 1 ,H''} = 8ttG 1 



dpi 



-8ttG 7 ' 



Pi 



(3.24) 
(3.25) 

(3.26) 



p[d cos(^ici) 


[sin(/2 2 c 2 ) - 


- sm(p' 3 c 3 )] 




2 7 A Pl 




jl' 2 C2 COs(p, 2 C 2 ) 


[sin(^ci) - 


- sin(/z' 3 c 3 )] 




2 7 Api 




p' 3 c 3 cos(p,' 3 c 3 ) 


[sin(/2ici) - 


- sm(p' 2 c 2 )] 


2 7 Api 



dpi 
dt" 



{p 1 ,Hp} = -8*G 7 d *^ 

p[ cos(/zici) [sin(/2 2 c 2 ) + sin(/2 3 c 3 )] 



(3.27) 



7A 



Note that substituting ( 2.51[ ) into (3.271 gives 
sm(p[ci) 



Mi 

7L1 [ H x 

2 \cos(/i 2 c 2 ) cos(/I 3 c 3 ) cos^ci) 



(3.28) 



Ho 



H\ + H2 Hi + H 3 



which is the modification of ( 2.53[ ) with quantum correc- 
tions. 



From ( |3.26| and (|3_27j), we have 
dc\ 



2 ^dt» 
which yields 



(>2 



dp2 

dt" 



- C3 



dP3 



167rG 7 p , (3.29) 



0, 



(3.30) 



— (p/c/ - pjCj) 
and hence, in accordance with the constant parameters 



used for classical solutions in (2.591, we set 
Pl ci = 8vrG 7 fi [AC/ + /(*)] 



(3.31) 



with the time- varying function fit). In the classical 
regime, pici become constant and so does fit). Taking 
(3.311 into the classical Hamiltonian constraint (2.581, 



we have 
*3 = 2((/C : 



-> f f{t)){IC 3 + f(t)) + cyclic terms J , (3.32) 



which gives the limiting values for fit) in the classical 
regime: 



/(*) 



or - -K 



as pjCj 



0. 



(3.33) 



That is, starting with f — (or / = — 2/C/3) in the clas- 
sical regime, / decreases (or increases) towards the big 
bounces and finally across the bounces / ends up with the 
other constant / = — 2/C/3 (or / = 0) on the other side 
of the classical regime. Equivalently, it can be said that 
the big bounces connect the classical solutions specified 
by K\, K2, k 3 on one asymptotic side with its "conjugate" 
counterpart on the other asymptotic side specified by the 
"conjugate" parameters: 



Kl,K 2 ,K3 



Kl 



K3 



(3.34) 



[Note th at the conj ugate parameters also satisfy the con- 
straints (I2.25I) and (12. 27b.] 



complicated constraint: 

'4(^2 



Substituting (|3.3l|) into (3.21 1, H'l, = then gives the 



2p 3 



f) 



p 3/2 

■ cyclic terms, 



4(^3 + /) 



,3/2 



where we define 



4 := SnG-fhA 1 / 2 



'pi 



(3.35) 



(3.36) 



for convenience. In particular, considering the classical 
limit: / -> or / -> -2 /C/3 and p 3 / 2 := y/pip 2 P3 > t\, 



we can show that (3.351 implies (2.261 



Substituting (3.311 into (3.27) and regarding <j> as the 



emergent time via (3.25), we have 
inGhp^ 2 



1 dpi 

Pi dcj> 



• 4 



p 



,3/2 



(3.37) 



4(^2 + /) 



P 



3/2 



4(^3 + /) 



p 



3/2 



Similar to the case of (|3.15|, p\ gets bounced once 
• )" term in 



"cos(- • • )" term in (3.37) flips sign. This happens when 

n3/2 



P 



3/2 







np 



(3.38) 



Assuming P2 and p 3 also get bounced roughly around the 
same moment, at which (3.381 is satisfied, we have the 
approximation: 



sin 



cos 



4(^2 + /) 



,3/2 



4(AC 2 -£i) 



,3/2 



1 - 



4(^2 -/Cl) 
p 3/2 

4(^2- /CQ 2 
2p 3 



(3.39) 
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and similar result for sin 

14(^3 + f)/p 3/2 )- Substitut- 
ing (3.38), (3.391 and (2.26) into ( 3.35[ ), we can solve: 



/>3 
7,3/2 



\K 1 -K 2 \ \Kx-K 3 \ 



X |^(k 2 + «2 + K-l) + (Kl - K 2 )(ni - K 3 ) 

-([( K l + K 2 + K l) + ( K l - « 2 )(ki - K 3 )]' 



-3(«i - k 2 ) 2 (ki - 



1/2" 



1/2 



=: K-iFfaKuKa). (3.40) 
The approximation is very accurate provided that 
F(ki;k 2 ,k 3 )\k2 - Kl| < 7T, 

i ? '(«i;K2i«3)|«3 -Kl| <7T 

4(/C 2 -/d) 2 4(/C 3 -/d) 2 



p3 



p3 



< IT 



(3.41) 
(3.42) 



The criterion (3.41) is satisfied for the "near-isotropic" 
solutions (i.e. ki, k 2 , K3 are not so different, or more pre- 



cisely, («i — n 2 y 



(ki-k 3 ) 2 < k 2 - 



-k\- 



j). In practice, 



however, we only need the left hand side of (3.41 ) to be 



fairly smaller than the right hand side in order to have 
(3.391 an good approximation, since the Taylor series of 
cos 2 = 1 — ir 2 /2 + • • • converges very rapidly. There- 
fore, for most solutions (even Kasner-like) , we can infer 
from (3.401 that that big bounces take place up to three 



times, once in each direction when the matter density p^ 
approaches one of the three critical values p/ iC rit given 
by the Planckian density ppj := (87rG7 2 A) _1 times a 
numerical factor: 



1 



Pl,crit ~ 7j K <p F («i;«2,«3)/9PJ- 



(3.43) 



For the solutions far from isotropic and thus (3.411 is 
badly violated, the above approximation is no longer 
good. Nevertheless, the big bounces still take place as 
most as three times at the moments when the "cos(- • • )" 
term flips sign in (3.37 ) . The matter density p^ is still the 



indication of the bounce occurrence, although the critical 



values could be quite different from (3.431. 

The differential equations ( 3.25 1— ( 3.27 1 can be solved 
numerically for given initial conditions! [See Footnote [7] 



for specifying the initial condition. Here, given pi(4>o), 
K-i and fCp, f((f>o) can b e obtained via (3.351 and then 
c/(0 o ) are fixed by (3.31 ).] The evolutions of pi (</)), a/(0) 



and P(/>(4>) are depicted in parts (a), (b) and (c) respec- 
tively in FIG.[3]for Kasner-unlikc solution and in FIG. [4] 
for Kasner-like solution. Note that in FIG.[3j the epochs 
of big bounces indicated by the dotted lines in (a) are 
very close to the moments corresponding to p/ jC rit indi- 
cated in (c). In FIG. |1J the moments indicated in (a) 
and those in (c) are quite close to eac h other but slightly 
different due to the fact that (3.411 is now only fairly 
satisfied. Furthermore, pi,p 2 ,p 3 bounce roughly around 
the same time, contrast to the p-scheme, in which the 
epochs of bounces in different directions can be arbitrar- 
ily separate. Also notice that each individual curve of 
Pi(<p) and ai(4>) in logarithmic scale becomes a straight 
line before and after the bounce, but the slope is changed, 
showing that the bounces conjoin two conjugate classical 
solutions. Finally, as in the p-schemc, the planar collapse 
in the Kasner-like case is not resolved. 

Contrast to the /I-scheme dynamics, in which the di- 
rectional density gj is the indication of bounces, it is the 
ordinary matter density that signals the occurrence 
of bounces in /i'-scheme. Unlike qi, the quantity of is 
independent of the choice of V since 



Pi 



V 2 



2pip 2 p 3 



V 2 



(3.44) 



Furthermore, (3.281 implies that the quantity p'jCi de- 
pends only on Hi, H 2 , H3 and is independent of V. Con- 



sequently, (3.31 1 and (3.35 1 tell us: For given initial phys- 
ical conditions <j>\t Q and Hj\ ta , K,<j,, Kj and f(t) all scale 
as oc V| to . Therefore, the effective dynamics given by 



(3.371 is completely independent of the choice of V as 
is the classical dynamics. In particular, the choice of V 
has no effect on the numerical factor k^F 2 (ki; kj, kk)/2 
appearing in (3.431. This is a desirable feature which 



/1-scheme does not have. [However, if we further impose 
the quantum corrections on the eigenvalue of the inverse 
triad operator as mentioned in (3.3 1, this invariance is 
broken again.] 

Even though /2'-scheme is independent of V, in case 
when V has a global meaning, the bounce condition (p^ f» 

Pi, crit 

) can be understood as: The physical volume of V 
(p 3 / 2 = V) gets bounced when it undergoes the Planck 
regime (times a numerical value 2k^ ) 2 F~ 2 (ki; kj, Kk)) 
measured by the reference of the momentum p^. 



IV. SCALING SYMMETRY AND RELATIONAL 
MEASUREMENTS 



With LQC discreteness corrections, both the /x- and p'- 
schemes show that the big bang singularity (both Kasner- 



like and Kasner-unlike) is resolved and replaced by the 
big bounces. The indication of the bounces is the direc- 
tional densities qi in /i-scheme whereas it is the matter 
density in p'-schcme. The detailed evolution of each 
individual pj shows that the bounces occur up to three 
times, once in each direction, when gj in p-scheme or p^ 
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(a) 



io 4 




1.0 1.5 2.0 2.5 3.0 
-(0-0„) <G-' 12 ) 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

-(0-0o) (C-" 2 ) 




FIG. 1: Kasner-unlike solution in /i-scheme effective dynamics, ki = 1/5, K2 = 1/4, K3 = 11/20, k^, — ^/ll9/200; 
pi((^o) = 1. x 10 5 4i, P2(0o) = 2. x 10 5 4i, Ps(0o) = 3. x 10 5 4i; and p 4> = 2.x IQphs/S^G (i.e., /Ck = 2. x 10 3 ). The red lines 



are for pi, 01, gi; green for j>2, 02, 52; and blue for ps, 03, £(3. The values of gj |Cr it given by (3.181 are pointed by the arrows 
in (c). (The Barbero-Immirzi parameter is set to 7 = 1.) 




lxlO 4 
5000 



0.5 1.0 1.5 

-(0-00) (G- |/2 ) 




0.5 1.0 1.5 

-(0-0o) (C-" 2 ) 



(c) 




0.5 1.0 1.5 

-(0-0,,) (G-" 2 ) 



FIG. 2: Kasner-like solution in /i-scheme effective dynamics, fti = —1/4, K2 = 3/4, K3 = 1/2, = 1/vo; pi(4>o) = 
1. x 10 5 4i, P2(^o) = 2. x 10 5 4i, Ps(0o) = 3. x 10 5 4b and = 2. x 10 3 /iv / 87rG (i.e., K,^ = 2. x 10 3 ). 



in //-scheme approaches its critical value. 

On the other hand, the planar collapse is not resolved 
but one of the length scale factors ai continues the van- 
ishing behavior in the Kasner-like case. This is ex- 



pected since the classical solutions (2.591 and (2.621 yield 
ftiCjjpfjCi — > (and also /j]cj — > in the /i -schcmc) 
toward the planar collapse and therefore the quantum 
corrections become more and more negligible. 

The fact that smallness of pi (not of ai) signals the oc- 
currence of big bounces seems to support the suggestion 
that "area is more fundamental than length in LQG" , 
although whether this is simply a technical artifact or 
reflects some deep physics is still not clear. (See Sec- 
tion VII. B of [16] for some comments on this aspect and 
[PT] for more details.) Meanwhile, as the length opera- 
tor has been shown to have a discrete spectrum [TH] , the 
fact that the vanishing of the length scale factor in the 
planar collapse is not stopped seems to contradict the 
discreteness of the length spectrum. Whether we miss 
some important ingredients when imposing the funda- 
mental discreteness of LQG in the LQC construction or 
indeed area is more essential than length remains an open 
question for further investigation. 

It has also been noted that the classical dynamics and 
the effective dynamics in /i'-scheme are both independent 
of the choice of the finite sized cell V, while the effective 



dynamics in /z-scheme depends on the physical size of V. 
This can be rephrased in terms of the scaling symmetry; 8 
that is, the classical dynamics and /i'-scheme effective 
dynamics are invariant under the following scaling: 



Pl,P2,P3 

ci, c 2 , c 3 
K.i 



hhpi, hkP2, hhP3, 

llCi, l 2 C 2 , I3C3, 

hhhPct» 



(4.1) 



[Note that the scaling for K-i should be accompanied by 
the same scaling on / in /i'-scheme; in this case, that 
is / — ► h^hf-] By contrast, the /i-scheme does not 
respect this scaling except for the special case with l\ — 
h = h- In this sense, /i-scheme has less symmetry of 
scaling than /i'-scheme and the classical dynamics. 



8 A dynamical system is said to be invariant under a certain scaling 
if for a given solution (pj(t), c/(i), rf>(t) and p^) to the dynamics, 
the rescaled functions also satisfy the equations of motion (i.e. 
Hamilton's equations and vanishing of Hamiltonian constraint). 
For cl assica l dyn a mics, the equations t o be s atisfie d are ( 2. 15 1 , 
\2.\G\ a nd j2.56H2.58t ; f or ^-s cheme, j3^- j3.1lt ; and kxpV- 
scheme, 13.241-13.271 and 13.35 1. 
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FIG. 3: Kasner-unlike solution in (j'-scheme effective dynamics. Ki = 1/5, K2 = 1/4, = 11/20, = ^/ll9/200; 
pi(^o) = 1. x 10 5 4i, P2(0o) = 2. x 10 5 4i, Pa{4>o) = 3. X 10 6 4i5 an <l = 2. x IO^a/SttG (i.e., /Ck,> = 2. x 10 3 ). The values of 
p/,crit given by (3.431 are pointed by the arrows in (c). The epochs of big bounces are indicated by dotted lines in (a), which 
are very close to the moments corresponding to p/, C rit, shown in (c). 




(c) 



0.4 0.6 0.8 1.0 
-(0-00) (C-" 2 ) 



0.4 0.6 0.8 
-(0-00) (C"" 2 ) 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 



FIG. 4: Kasner-like solution in //-scheme effective dynamics, ki = —1/4, /t2 = 3/4, K3 = 1/2, k$ — l/v8; pi(4>o) = 
1. x 10 4 4i, MM = 2. x 10 4 4i, P3(<M = 3. x 10 4 4b and P4> = 2. x lO^v 7 ^? (i.e., = 2. x 10 3 ). The epochs of big 
bounces are only slightly different from the moments corresponding to pi jC rit. 



Apart from the scaling related to the dependence on 
V, all three theories (classical, /I-scheme, /i'-scheme) also 
admit the scaling symmetry under: 



Pi,P2,Pa 

Cl,C 2 ,C 3 

P<1> 



\<f>, 

X 2/3 P1, X 2/Z P2, X 2/3 P3, 

AV3 Cl , X^c,, AV3 C3) 

XlCfj,, 
XKi. 



(4.2) 



This is reminiscent of the idea as suggested in [HO 
20 that length/ area/ volume is measurable only if the 
line/surface/bulk is coupled with the material reference. 

Another symmetry of scaling respected by the classical 
dynamics is given by 



t ► T)t, 

7 — ► £7, 

ci,c 2 ,c 3 — ► £,n~~ 1 c 1 , £r)~ 1 c 2) £J7 _1 C 3 , 

— ► »7 _1 |ty> 

^0 — ► V~ fc<t>i 

Ki — ► r\ Ki. 



The scaling symmetry regarding 7 — > £7 is expected, 
since the Barbero-Immirzi parameter 7 has no effect on 
the classical dynamics. The scaling symmetry regarding 
t — ► rjt is also easy to understand, since there is no 
temporal scale introduced in the Hamiltonian. However, 
very surprisingly, the scaling symmetry involving t — > 
■qt is violated for both /2-scheme and /i'-scheme effective 
dynamics. Curiously, this symmetry is restored if t — ► 
•qt is accompanied by 7 — ► £7 and one extra scaling is 
also imposed at the same time: 



r 4 v 



2 - 2 A. 



(4.4) 



(4.3) 



This intriguing observation seems to suggest, albeit spec- 
ulatively, that in the context of quantum gravity the fun- 
damental scale (area gap) in spatial geometry gives rise to 
a temporal scale via the non-locality of quantum gravity 
(i.e., using holonomies) and the Barbero-Immirzi param- 
eter 7 somehow plays the role bridging the scalings in 
time and space. [This reminds us that, in LQG, the pre- 
cise value of the area gap A is proportional to 7, and 7 
is also the parameter which relates the intrinsic geome- 
try (encoded by spin connection T l a ) with the extrinsic 
curvature (K a l ) via Aj = T i a — 7 AY .] 

Meanwhile, related to the above observations, the 
physical meaning of the directional factors qi and mat- 
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ter density p^ can be interpreted as the (inverse of) area 
and volume scales, measured by the reference of the mat- 
ter content. In this regard, we may say that the big 
bounces take place when one of the area scales or the 
volume scale becomes very small by the reference of the 
matter momentum. It is then tempting to regard not 
only (j) as the "internal clock" (emergent time) but also 
p<f, as the "internal rod" — namely, the measurement of 
both temporal and spatial geometries makes sense only 
in the presence of matter content. 

The observations about the symmetry of the scalings 
(j) — ► \<j), t — ► -qt and 7 — > £7 may support the ideas of 
the relational interpretation of quantum mechanics with 
real rods and clocks such as studied in [21] (see also [19j 
20 ), although the link is far from clear. However, as 
the theories studied in this paper are highly simplified, 
these observations could be just artifacts and we should 
not take this speculation too seriously until this aspect 
is further investigated for more sophisticated models. 



Unfortunately, all the scaling symmetries break down 
in the detailed construction of LQC with /t-scheme (the 
strategy to construct the full theory of LQC in //'-scheme 
is still not clear) even for the isotropic model (where p- 
and /i'-schcmcs arc identical). The full quantum theory 
only respects the scaling symmetries at the leading order. 
This is due to the fact that the quantum evolution in the 
full theory of LQC is governed by a difference equation, 
in which the step size of difference introduces an addi- 
tional scale in the deep Planck regime (see [7] for Bianchi 
I model and [I] for the isotropic model). In fact, already 
in the level of effective dynamics, the scaling symmetries 
are violated if we further take into account the LQC cor- 
rections on the inverse triad as given by (3.3 1. For the 



full theory of LQC, if we take the aforementioned sym- 
metries seriously, we might be able to revise the detailed 
construction in the spirit of relational quantum theory 
such that the step size in the difference equation scales 
adaptively by the reference of the matter content. 



V. SUMMARY AND OUTLOOK 

To summarize, we list the important facts for the 
classical dynamics, /i-scheme effective dynamics and p!- 
scheme effective dynamics in TABLE [I] In the following, 
the main results are restated and some feasible extensions 
are remarked. 

With the LQC discreteness corrections, the effective 
dynamics shows that the classical singularities (both 
Kasner-like and Kasner-unlike) are resolved and replaced 
by the big bounces, which take place up to three times, 
once in each diagonal direction. In /2-scheme, it is the 
directional densities qj that signal the occurrence of big 
bounces whenever each of qi approaches its critical value. 
In /i'-scheme, the indication is the matter density p^ and 
the big bounces happen at the moments when p$ ap- 
proaches one of the three critical values. 

Furthermore, the detailed evolution shows that the 
equations of motion (in terms of emergent time <f>) for pi 
in different diagonal directions are completely decoupled 
and evolve independently of one another in /I-scheme; as 
a result the moments of three bounces can be arbitrarily 
separate. By contrast, in /I'-scheme, equations of mo- 
tion for pi are coupled through p^ and thus p\, P2, and 
P3 bounce roughly around the same moment. In both 
schemes, there are three bounces in general unless the 
initial conditions are delicately fine tuned. 

On the other hand, the planar collapse is not resolved 
but one of the length scale factors ai continues the van- 
ishing behavior in the Kasner-like case. Whether this 
suggests that area is more fundamental than length or 
conflicts with the discrete spectrum of length operator in 
LQG requires further investigation. 

Across the bounces, the equations of motion again 
come closer and closer to the classical ones. Hence, the 



semiclassicality is retained on both asymptotic sides of 
the evolution. Given a classical solution in one asymp- 
totic side, the evolution ends up with another classical 
solution on the other asymptotic side, which is either an- 
tipodal or conjugate to the given solution in p / p' -scheme 
respectively. 

In regard to the finite sized cell V chosen to make sense 
of the Hamiltonian formalism, the effective dynamics in 
/2-scheme depends on the choice of V, and thereby re- 
acts to the macroscopic scales introduced by the bound- 
ary conditions. (In terms of symmetry, it is said that 
/I-scheme has less scaling symmetry than /2'-scheme and 
classical dynamics.) The effective dynamics in /t'-scheme, 
by contrast, is completely independent of V as is the clas- 
sical dynamics. (The issue of dependence on V may be 
related to the instability of /2-scheme indicated in [T5].') 
In case that the physical size of V has a global mean- 
ing (such as in the compactified Bianchi I model or the 
finite sized homogeneous patches in BKL scenario), the 
condition for the bounce occurrence can be rephrased: In 
p,/ pi '-scheme (respectively), the physical area/volume of 
the surfaces/bulk of V gets bounced when it undergoes 
the Planck regime (times a numerical value) measured by 
the reference of the momentum 

While the /t'-scheme has the advantage that its effec- 
tive dynamics is independent of V, the full theory of LQC 
based on /2'-scheme is difficult to construct. Both p- and 
/t'-schemes have desirable merits and it is still disputable 
which one is more faithful implementation of the under- 
lying physics of quantum geometry. 

In addition to the symmetry related to the choice of 
V, both schemes admit additional symmetries of scal- 
ing, which are reminiscent of the relational interpretation 
of quantum mechanics, featuring the ideas of real rods 
and clocks. Furthermore, the symmetry involving the 
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classical dynamics 


effective dynamics in /i-scheme 


effective dynamics in //-scheme 


Pep = tly/ OTYLrl\,^ = V cp 


Pep = raVO7T(_r/L0 = \ <p 


P(p = ftVO7T(_r/L0 = V Cp 


pi = L 2 L s (a 2 a 3 ) = L 2 L 3 


pi = L 2 L 3 (a 2 a3) = L 2 L 3 


pi = L 2 L 3 {a 2 a 3 ) = L 2 L 3 


cj = yLiHi = 7L1 




sinf^ci) 

Pi _ 

7L1 ( H 1 +H 3 | H 1 +H 2 H 2 +H 3 \ 


2 \ cos(/^2C2) cos(^3C3) cos(/^ici)J 


2 cos(^2 c 2) Cos(/LigC3) COs(MiCx) J 


pici = 8nGyhlCi = ^WHi 


Mi 


p/ C/ = 87^7/1 [/C/ + /(*)] 
/ — ► 0, — in classical regime 


/C| = 2 (IC2IC3 + IC1K.3 + K.itC 2 ) 


Kj> = 2 (JC2JC3 + /C1/C3 + /C1/C2) 


/C| = 2 (K.2IC3 + AC1AC3 + /C1/C2) and 
^ = 2{^sin(^)sin(^) 

+ cyclic terms | 


1 d Ji = 87rGft (^2 + /c 3 ) 


pT^-^ COS ^ Ci )^+^) 


1 dpi _ SttGH co „ 
pi dtp p,f, " y p 3/2 y 

x«£ {sin(^™') + sin(^^)} 


pi -> 

toward singularity 
at the same moment 


bounces whenever 

2 2 

Qi ■= jps = 7IPP1; 
epochs of the bounces can 
be arbitrarily separate 


pi bounces at the moment when 

2 

: = apiptpa ~ 2^ («i;«2,k3)ppi; 

pi,P2,P3 bounce roughly around 
the same moment 


no big bounce; 
ki fixed 


the big bounces conjoin the pair 
of "antipodal" classical solutions: 

K[ < > — Kl 


the big bounces conjoin the pair 
of "conjugate" classical solutions: 

Kl < > Kl - | 


symmetry of scaling: 

* > Tjt 

7 — > £7 
</> — > A</> 

Pi,P2,P3 — >hh\ 2/3 pi, ■■■ 
ci,c 2 ,c 3 — > ZiA 1/3 £?7 -1 ci, ... 
P0 — ► hhhXr)~ 1 P4, 
— > lihh^il 1 lC</> 
KLi — > hhhXr]^ 1 ICi 


symmetry of scaling: 

* > Tjt 

7 — ► £7 

4> — ► \<t> 

Pl,P2,P3 ► l 2 X 2/3 pi, . . . 

Cl,C 2 ,C 3 ► l\ 1 ' 3 ^ 1 Cl, ... 

p<p — ► l 3 \ri~ 1 p r p 
K.4 — > l''\ri~ 1 K, r p 

id — > fx-n^tCi 

A — > £- 2 rf& 


symmetry of scaling: 

i — > r?t 

7 — ► £7 
4> — > \</> 

Pi,p 2 ,P3 — > khX 2/3 Pi, ••• 
ci,c 2 ,c 3 — > iiA 1/3 £r?~ 1 ci J ... 
— > hhhXri^pj, 

K<p > lil 2 ls\r]~ 1 IC,p 

(Ki + /) — ► hhkXr]- 1 ^! + /) 
A-^C" 2 ^ 2 A 



TABLE I: Summary of the classical dynamics, /i-scheme effective dynamics and //-scheme effective dynamics. 



Barbero-Immirzi parameter is suggestive that the funda- 
mental scale (area gap) in spatial geometry may give rise 
to a fundamental scale in temporal measurement. These 
symmetries however break down in the construction for 
the full theory of LQC. 

A few possible extensions following the treatment of 
this paper seem to be within reach. First, as this paper 
focuses specifically on the model with a massless scalar 
field, it should be straightforward (with necessary ap- 
proximation) to extend the results to the models with 
inclusion of generic matters, which will fulfill the inves- 
tigation in |10j for the //-scheme. Another extension in- 
volves the inclusion of nontrivial potentials for the scalar 
field. In both extensions, most of the observations we 
obtained here are to be expected and the investigations 
in broader context would further support or oppose some 
of our speculations. 

Second, it would be instructive to compare the re- 



sults of this paper with the perturbative treatment of 
anisotropies, since it has been suggested that the big bang 
singularity is not resolved if anisotropies are treated as 
perturbations of an isotropic background [22] . Incorpo- 
ration of inhomogeneities, on the other hand, will be sub- 
stantially more difficult. In order to determine whether 
jl- or //-scheme is more physically sensible, we have to in- 
corporate inhomogeneities at least approximately or per- 
turbatively in the background of Bianchi I type. The 
perturbative treatment such as developed in or the 
strong gravity approximation proposed in |24j might be 
useful for this purpose. 

Finally, the results in principle could be extended to 
the Bianchi IX model. Since Bianchi IX model, along 
with Bianchi I model, plays a pivotal role for the infinite 
fragmentation of the homogeneous patches close to the 
singularity in the BKL scenario, the analysis may provide 
hints to how the loop quantum effects resolve the classical 
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singularities in generic situations and in return it may 
help arbitrate the tension between p- and /i'-schemes for 
the Bianchi I model. 
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dci 
~dF 



dpi 
dt' 



d H 1 

_ x sinQuJci) 



Mi 

M2 



sm(/i|c 3 ) 

P2 H f P3 

^3 



d H' 



dci 



7 pi cos(/x°ci) 

X 



sin(^c 2 ) sin(/^c 3 ) 

Vi H t, P3 



f-2 



f>3 



(A4) 



(A5) 



APPENDIX A: EFFECTIVE DYNAMICS IN 
^-SCHEME 

One of the virtues of the improved strategy (/I-scheme) 
in the isotropic model is to fix the serious drawback in the 
old precursor strategy (/i -scheme) that the critical value 
of the matter density at which the bounce occurs can be 
made arbitrarily small by increasing the momentum p^, 
thereby leading to wrong semiclassical behavior. 

In the Bianchi I model, either the directional densities 
Qi (in /2-schcmc) or the matter density p^ (in /i'-schcmc) 
plays the same role as p^ does in the isotropic case. Hav- 
ing learned from the isotropic case, we expect that the 
critical values of qi or p^ at which the bounces occur 
can be made arbitrarily small by increasing the momen- 
tum p^ in /io-scheme but is independent of p^ in p- or 
/i'-scheme. The latter is what has been shown in the 
main text of this paper. 9 For comparison, the effective 
dynamics in /i -scheme is presented here. 10 

In the effective theory of /^ -scheme, we take the pre- 
scription to replace ci by sin(/if c/)//if with the fixed 
numbers p°j for discreteness. Analogous to (3.4), we have 
the effective (rescaled) Hamiltonian constraint: 



TTl _ r<p 



(Al) 



87rG 7 2 



sin(^c 2 ) sin(//gc 3 ) 



P2P3 + cyclic terms 



Again, the equations of motion are given by the Hamil- 
ton's equations and the constraint that the Hamiltonian 
must vanish: 



dp<f, 

IF 
d4> 

~dl' 



= {p^H , fl J = 



Pcf, is constant, (A2) 
(A3) 



9 The critical values £/ iCr it depend on p^ only through the ra- 
tios in /i-scheme; p/ iCr it depends on p^ only through 
k^F 2 (kj; kj, kk)/2 in /i'-scheme. 
10 For the Bianchi I LQC, the effective dynamics in /t -scheme was 
first studied in [8] for the vacuum solution. 



as well as 



^ o (c/,Pz) = => 

1 f sin(/4!c2) sin(^3C 3 ) 



4irG~/ 2 



^3 



(A6) 

P2P3 + cychc terms 



From ( A4 1 and ( A5 ) , we have 
shi(/^cj) 



d 
df 



Pr 



I'i 



which gives 



sin(^c 7 ) 

pi = SuG^hlCi. 

M? 



(A7) 



(A8) 



Taking ( A8 ) into ( A6 I gives the same constraints on the 



constant parameters as in (2.261 



Substituting (A8) into (A5l yields 



1 dpi 
pi dt' 



8-kGH cos(^ci)(/C2 + £3)- 



By (A3) and cos a; = ±\/T — sin 2 x, (A9) leads to 

/ \ 2 / 3 



11 = ^^-' 

Pi dcj) 



1 



Qi 



Ho 

@L crit 



(A9) 



1/2 



(A10) 

which gives the bouncing solutions with the behaviors 
similar to those given by (3.161 except that the critical 



value of qi at which the big bounce takes place is given 

by 



Mo 
Si, crit 



Kl 



-1 3/2 



1 



(All) 



which can be made arbitrarily small by increasing the 
value of Ptj, . As a result , /i -scheme gives wrong semiclas- 
sical behavior and should be improved by p- or /i'-scheme 
to fix this problem. 
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APPENDIX B: QUANTIZATION IN THE 
/i-SCHEMES 

In this appendix, we give a heuristic argument for the 
prescription given by (3.11, starting from the Hamilto- 



nian constraint of LQG. The motivations for both p,- and 
/x'-schemes are addressed in detail. The advantages and 
drawbacks of both schemes are also remarked. 

The gravitational part of the classical Hamiltonian 
constraint in the full theory of LQG is given by 



H,. 



grav 



' 1 <PxNe- 1 \e i > k F* b &' J E? k 



8ttG 



-2{\ + 1 2 )K [a l K b] ] E a l E b 3 ^ 
= J d 3 xNe-W k F* b E a j E b k , (Bl) 

where e := | detE] 1 / 2 and in the last line we exploit the 
fact that for homogeneous and spatially flat models the 
two terms inside the curly bracket are proportional to 
each other (because —jK a l — A a l as the spatial slice is 
fiat and F* b — e 1 j k A a 3 A^ 3 as it is homogeneous). Fur- 
thermore, the lapse N can be assumed to be constant 
because of homogeneity and we set N = 1 . 

When applied to Bianchi I models, the integral is re- 
stricted to a finite sized cell V as prescribed by (2.6 1. 



Accordingly, we should also adopt the replacement rules 
(2.391 and (2.401. This procedure gives 



1 



8ttG7 2 



x ^ k e\ m A p l A q m EP J E\) 



8ttG7 2 



and 



x (C2C3P2P3 + cyclic terms) 
(C2C3P2P3 + cyclic terms) 
8ttG 7 2 e V 



|det£| = V- 3 \det(E j l °e c ;)\ 



(B2) 



V~ 



P1P2P3 



y 2 \piP2P3\ 



Consequently, we have 



(C2P2C3P3 + C1P1C3P3 + C 1 p 1 C 2 P2) 

8ttG7 2 v /|piP2P3| 



(B3) 



(B4) 



which is exactly the same as H gT3lV given in (2.54) 



When the quantization is performed in the context of 
LQC, there are two loop quantum corrections. The first is 

the modificat ion on the inverse triad operator as 
indicated in (3.3 1, which is negligible and ignored in this 



Aa (or a) do not exist and should be replaced with 
holonomies (or exponentials of c/). 

Following the standard techniques in gauge theories, 
the curvature component F* b can be expressed in terms 
of holonomies (i.e. Wilson loops). Given a small surface 
a center in x, Stokes' theorem allows us to write the 
curvature component as 



nFUx) 



1 



c ab Jot 



TiFi d dx c A dx d 



V exp 



da 



(B5) 



where 2iti = oi are the Pauli matrices; da is the bound- 
ary loop of a; and the coordinate area of a projected in 
"a6-direction" is given by 



eZ b = I dx a Adx b . 



(B6) 



This is a good approximation provided e" b is small 
enough and in fact it becomes exact in the continuous 
limit e a ab - 0. 

For the Bianchi I model, we choose a = Djk to be a 
small rectangular surface parallel to the J-K plane. The 
coordinate lengthes of the edges of a are denoted as JIjLk 
and P>kLk and accordingly the coordinate area is given 
by e°j K = = JxjJikLjL k (see FIG. [H). Applying 
(2.401 and noting that F l ab = e i jk A a l A b j gives 



F, 



ah 



Fjk '■— tIJK(LjL k ) 1 CjCk, 



we read off from ( B5 1 that 



Tr 



Z JK 



Tl 



\ n JK 



(B7) 



(B8) 



where (I,J,K) is any permutation of (1,2,3) and the 
holonomy along the edges of Djk 1S: 



(B9) 



with lif" 1 ^ being the holonomy along the individual edge: 



hf l] := Vexp 



1 ^ 1 



T iA a dx 



exp (/i/C/r/) 



cos[ ^ r)+2s in [ — ) r I , 



(BIO) 



which follows 
Tr \ Tl h^ J ' fiK) 



(-UK 



sm(fijcj)sin(fi K c K )- (Bll) 



Putting ( |B7| |, ( p8| and ( |B11| altogether, we then have 
sm(fljcj) sm(fi KCK ± 



cjc K 



paper. The second is due to the fact that the connections which is exactly the prescription used in (3.1 ). 
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If Ojk shrinks to a point, this gives the continuous 
limit /ij — > and we recover the classical Hamiltonian 



constraint in (B4). However, the very feature of LGC is 



that the continuous limit does not exist and the failure of 
the limit to exist is intimately related with the underly- 
ing quantum geometry of LQG, where eigenvalues of the 
area operator are discrete. In LQC, to implement the dis- 
creteness as imprint from the full theory of LQG, we have 
to set jli to be finite. The question now is: What values 
should fli be set to faithfully reflect the fundamental dis- 
creteness of quantum geometry? Two strategies (/2- and 
//-schemes) are discussed below. 

Since the eigenvalue spectrum of the area operator has 
an area gap A, to impose the discreteness, it is straight- 
forward to set the physical area of □ jk to be A (depicted 
in FIG.[5ja); i.e. 



a J a K e JK — a J a K^J^KLjL K 

= \tI.JKPl\ fij^K = A, 



(B13) 



which gives the "//-scheme" in (3.6 1 



However, a closer consideration seems to suggest that 
equating the physical area of Ojk to A might not be 
completely sensible. Recall that the validity of ( |B5| is 
based on Stokes' theorem, which does not invoke the 
metric of the surface a at all. The smallness of a to 



ensure good approximation of ( B5 1 does not directly re- 
fer to metric smallness, either. Furthermore, if we think 
of the picture of spin networks in LQG, what associated 
with the physical areas are links (holonomics) of the spin 
network; the area enclosed by the loop of holonomy links 
is irrelevant! Therefore, back to LQC, instead of shrink- 
ing the area of Ojk, we should associate each edge of 
Ojk with an area and then shrink the associated areas 
to A. (See FIG.[5]b.) The edge of O jk in ./-direction is 
of coordinate length p,jLj, with which, most naturally, 
we associate a rectangle Oj parallel to I-K plane of co- 
ordinate lengths fijLi and JIjLk on its edges. (That is, 
the ratio of the edge of Oj to the edge of V in the //in- 
direction is set to be the same as the ratio of the edge of 
Ojk to the edge of V in the ./-direction.) We then set 
the physical area of Oj to A; i.e. 



ajaK^jLjhK = l-i 2 j\pj 



(B14) 



which gives the "/i-scheme" in (3.51. 

[Note that in both pictures or p'- and /i-schemes, if we 
rather set the coordinate area of Ojk or Oj to be A, 
we end up with constant /i/. This is the "/i -schcmc" 
used in the precursor strategy. The wrong semiclassical 
behavior of the /x -scheme originates from the problem 
that the coordinate area (both of Ojk and of Oj) is not 
totally physically relevant.] 

In the full theory of LQG, when the Hamiltonian acts 
on an spin network state, it adds a new link of spin-1/2 
and the coloring of the links on which the new link is at- 
tached is increased or decreased by 1/2 (see FIG.[6]a-b). 
The resulting spin network state is equivalent to the orig- 
inal state superimposed with a triangular loop of spin-1 /2 



links (FIG.|6]c). This triangular loop is essentially the 
Wilson loop discussed above. From this perspective, the 
//-scheme can be understood as associating the area in 
pink in FIG.[6jd with A, whereas the /i-scheme as associ- 
ating the area in blue in FIG. |6]e with A. In the context 
of spin network states, the coloring of a link corresponds 
to the area of the surface to which the link penetrates and 
the smallest coloring spin-1/2 gives rise to the area gap 
A. Therefore, it is in this sense that /i-scheme is a more 
direct implementation of the underlying discreteness of 
quantum geometry than //-scheme. 

Furthermore, in the full quantum theory of LQC, /2- 
scheme has the important virtue that we can define the 
affiinc variables vj via 



d 
dvj 



Pi 



_d_ 
dpi 



(B15) 



such that Hamiltonian constraint of the full quantum the- 
ory gives the evolution as a difference equation in terms 
of vi and therefore the methodology used for the isotropic 
LQC can be easily applied [7]. This strategy fails in p,'- 
scheme since 



_d_ 
dpi 



/'./ 



_d_ 
dpj 



^0 



(B16) 



and hence the corresponding affine variables do not ex- 
ist. This makes it difficult to construct the full quantum 
theory of Bianchi I LQC in //-scheme. 

On the other hand, as studied in Section III //-scheme 



has the advantage over /i-scheme that the effective dy- 
namics in //-scheme is independent of the choice of V. 
The difference for this point between p- and //-schemes 
can be understood, heuristically but instructively, by es- 
timating the quantity /i/cj with the classical formulae; 
that is, plugging (2.531 for a, we have 



i lCl « 7 A 1/2 #i ( S 



L2L3 



1/2 



/2'jc/ « 1 & 1 ' 2 Hi. 



(B17) 
(B18) 



Since the quantity /i/c/ indicates how significant the 
quantum corrections are (quantum corrections are neg- 
ligible if piici <C 1 ), (B17) tells that in /i-scheme, the 



place at which the quantum effects become appreciable 
is tied up with not only the Hubble rates Hi but also the 
physical geometry of V. By contrast, (B18l shows that 



the different choice of V is irrelevant in //-scheme. 

In the language of [TS] , /i-scheme corresponds to a lat- 
tice refinement model whereby the number of vertices is 
proportional to the transverse area. A stability analysis 
suggests that /i-scheme leads to an unstable difference 
equation in the full theory of LQC and thus may not 
represent a good quantization scheme. 

Both /i- and //-schemes have desirable and undesirable 
features of their own. In order to understood them more 
deeply, in the main text we study both schemes and their 
ramifications at the level of effective dynamics. 
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FIG. 6: (a)-(b) The action of Hamiltonian acting on a spin network state, (c) The resulting state is equivalent to the original 
one superposed with a triangular loop, (d) The idea of p'-scheme. (e) The idea of /i-scheme. 
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